Key mRNAs and lncRNAs of pituitary that affect the reproduction of FecB + + small tail han sheep.

BACKGROUND: The pituitary directly regulates the reproductive process through follicle-stimulating hormone (FSH) and luteinizing hormone (LH). Transcriptomic research on the pituitaries of ewes with different FecB (fecundity Booroola) genotypes has shown that some key genes and lncRNAs play an important role in pituitary function and sheep fecundity. Our previous study found that ewes with FecB + + genotypes (without FecB mutation) still had individuals with more than one offspring per birth. It is hoped to analyze this phenomenon from the perspective of the pituitary transcriptome. RESULTS: The 12 Small Tail Han Sheep were equally divided into polytocous sheep in the follicular phase (PF), polytocous sheep in the luteal phase (PL), monotocous sheep in the follicular phase (MF), and monotocous sheep in the luteal phase (ML). Pituitary tissues were collected after estrus synchronous treatment for transcriptomic analysis. A total of 384 differentially expressed genes (DEGs) (182 in PF vs. MF and 202 in PL vs. ML) and 844 differentially expressed lncRNAs (DELs) (427 in PF vs. MF and 417 in PL vs. ML) were obtained from the polytocous-monotocous comparison groups in the two phases. Functional enrichment analysis showed that the DEGs in the two phases were enriched in signaling pathways known to play an important role in sheep fecundity, such as calcium ion binding and cAMP signaling pathways. A total of 1322 target relationship pairs (551 pairs in PF vs. MF and 771 pairs in PL vs. ML) were obtained for the target genes prediction of DELs, of which 29 DEL-DEG target relationship pairs (nine pairs in PF vs. MF and twenty pairs in PL vs. ML). In addition, the competing endogenous RNA (ceRNA) networks were constructed to explore the regulatory relationships of DEGs, and some important regulatory relationship pairs were obtained. CONCLUSION: According to the analysis results, we hypothesized that the pituitary first receives steroid hormone signals from the ovary and uterus and that VAV3 (Vav Guanine Nucleotide Exchange Factor 3), GABRG1 (Gamma-Aminobutyric Acid A Receptor, Gamma 1), and FNDC1 (Fibronectin Type III Domain Containing 1) played an important role in this process. Subsequently, the reproductive process was regulated by gonadotropins, and IGFBP1 (Insulin-like Growth Factor Binding Protein 1) was directly involved in this process, ultimately affecting litter size. In addition, TGIF1 (Transforming Growth Factor-Beta-Induced Factor 1) and TMEFF2 (Transmembrane Protein With EGF Like And Two Follistatin Like Domains 2) compensated for the effect of the FecB mutation and function by acting on TGF-ß/SMAD signaling pathway, an important pathway for sheep reproduction. These results provided a reference for understanding the mechanism of multiple births in Small Tail Han Sheep without FecB mutation.

Screening of Litter-Size-Associated SNPs in NOX4, PDE11A and GHR Genes of Sheep.

In previous studies, NOX4, PDE11A and GHR genes have been screened as important candidate genes for litter size in sheep by using the GWAS method; however, neither their effects on litter size nor the loci associated with litter size have been identified. In this study, three candidate loci (c.1057-4C > T in NOX4, c.1983C > T in PDE11A and c.1618C > T in GHR) were first screened based on our previous resequencing data of 10 sheep breeds. After the three loci were genotyped using Sequenom MassARRAY technology, we carried out population genetics analysis on the three loci and performed association analysis between the polymorphism of the three loci and the litter size of sheep. The results of population genetics analysis suggested that c.1057-4C > T in NOX4 and c.1983C > T in PDE11A may be subject to natural or artificial selection. The results of association analysis indicated that litter size was significantly associated with c.1057-4C > T in NOX4 and c.1983C > T in PDE11A (p < 0.05) in Small Tail Han sheep, and there was no significant interaction effect between the two loci on the litter size. In summary, c.1057-4C > T in NOX4 and c.1983 C > T in PDE11A can be considered candidate molecular markers for improving litter size in sheep.

Thyroid transcriptomic profiling reveals the differential regulation of lncRNA and mRNA related to prolificacy in small tail han sheep with FecB BB genotype.

BACKGROUND: The thyroid gland is an important endocrine gland in animals that secretes thyroid hormones and acts on various organs throughout the body. lncRNAs are long non-coding RNAs that play an important role in animal reproduction; however, there is a lack of understanding of their expression patterns and potential roles in the thyroid gland of the Small Tail Han (STH) sheep. In this study, we used RNA-Seq technology to examine the transcriptome expression pattern of the thyroid from the luteal phase (LP) and follicular phase (FP) of FecB BB (MM) STH sheep. RESULTS: We identified a total of 122 and 1287 differential expression lncRNAs (DELs) and differential expression mRNAs (DEGs), respectively, which were significantly differentially expressed. These DELs target genes and DEGs can be enriched in several signalling pathways related to the animal reproduction process. CONCLUSIONS: The expression profiles of DELs and DEGs in thyroid glands provide a more comprehensive resource for elucidating the reproductive regulatory mechanisms of STH sheep.

Transcriptomic Analysis Reveals Differentially Expressed Circular RNAs Associated with Fecundity in the Sheep Hypothalamus with Different FecB Genotypes.

Circular RNAs (circRNAs) are a specific type of noncoding RNA, and some have defined roles in cellular and biological processes. However, little is known about the role of circRNAs in follicular development in sheep with FecB (fecundity Booroola) mutations. Here, the expression profiles of circRNAs were investigated using RNA sequencing (RNA-seq) in the follicular phase (F) and the luteal phase (L) of FecB mutant homozygous (BB) and wild-type (WW) Small Tail Han sheep. A total of 38,979 circRNAs were identified, and 314, 343, 336, and 296 of them were differentially expressed (DE) between BB_F and BB_L, WW_F and WW_L, BB_F and WW_F, and BB_L and WW_L, respectively. The length, type, and chromosome distribution of the circRNAs and the expression characteristic between the circRNAs and their host genes in the sheep hypothalamus were ascertained. Enrichment analysis showed that the host genes of DE circRNAs in the follicular and luteal phases were annotated to MAPK, gap junctions, progesterone-mediated oocyte maturation, oocyte meiosis, and other hormone-related signaling pathways, and the different FecB genotypes were annotated to the gap junctions, circadian entrainment, MAPK, and other hormone-related signaling pathways. The competing endogenous RNA network prediction revealed that the 129 target miRNAs might be bound to 336 DE circRNAs. oar_circ_0000523 and oar_circ_0028984, which were specifically expressed during the follicular phase in the BB genotype sheep, probably acted as miRNA sponges involved in the regulation of LH synthesis and secretion. This study reveals the expression profiles and characterization of circRNAs at two phases of follicular development considering different FecB genotypes, thereby providing an improved understanding of the roles of circRNAs in the sheep hypothalamus and their involvement in follicular development and ovulation.

Thyroid transcriptomic profiling reveals the differential regulation of lncRNA and mRNA related to prolificacy in Small Tail Han sheep with FecB++ genotype.

The thyroid gland is an important endocrine gland in animals, which mainly secretes thyroid hormones and acts on various organs of the body. Long-chain non-coding RNA (lncRNA) plays an important role in animal reproduction. However, there is still a lack of understanding of their expression patterns and potential roles in the thyroid of Small Tail Han (STH) sheep. In this study, RNA-seq was used to examine the transcriptome expression patterns of lncRNAs and mRNAs in the follicular phase (ww_FT) and luteal phase (ww_LT) in FecB++ genotype STH Sheep. A total of 17,217 lncRNAs and 39,112 mRNAs were identified including 96 differentially expressed lncRNAs (DELs) and 1054 differentially expressed mRNAs (DEGs). Functional analysis of genes with significant differences in expression level showed that these genes could be enriched in Ras signalling pathway, hedgehog (HH) signalling pathway, ATP-binding cassette (ABC) transporters and other signalling pathways related to animal reproduction. In addition, through correlation analysis for lncRNA-mRNA co-expression and network construction, we found that LNC_009115 and LNC_005796 trans target NIK-related kinase (NRK) and poly(A)-specific ribonuclease (PARN). LNC_007189 and LNC_002045 trans target progesterone-induced blocking factor 1 (PIBF1), LNC_009013 trans targets small mothers against decapentaplegic (SMAD1) are related to animal reproduction. These genes add new resources for elucidating the regulatory mechanisms of reproduction in sheep with different reproductive cycles of the FecB++ genotype STH sheep.

Association Analyses between Single Nucleotide Polymorphisms in ZFAT, FBN1, FAM184B Genes and Litter Size of Xinggao Mutton Sheep.

Previous studies have screened key candidate genes for litter size in sheep, including fibrillin-1 (FBN1), family with sequence similarity 184 member B (FAM184B) and zinc finger and AT-hook domain containing (ZFAT). Therefore, it is necessary to verify these genes in the Xinggao mutton sheep population and determine the associated loci for litter size. In this study, three loci (FBN1 g.160338382 T > C, FAM184B g.398531673 C > T and ZFAT g.20150315 C > T) were firstly screened based on the population differentiation coefficient between the polytocous and monotocous sheep groups. Then, population genetic analysis and association analysis were performed on these loci. The results revealed that the g.160338382 T > C in FBN1 was significantly associated with the litter size of sheep. Moreover, there was no significant interaction effect between the g.160338382 T > C locus and FecB on litter size. Notably, g.160338382 T > C is adjacent to the anterior border of exon 58 and belongs to a splice polypyrimidine tract variant, which may lead to alternative splicing and ultimately cause changes in the structure and function of the protein. In summary, our results provided a potentially effective genetic marker for improving the litter size of sheep.

Thyroid Transcriptomics Revealed the Reproductive Regulation of miRNA in the Follicular and Luteal Phases in Small-Tail Han Sheep with Different FecB Genotypes.

MicroRNA (miRNA) is a type of endogenous short-stranded ncRNA that influences many biological processes such as animal growth, development and metabolism. The thyroid gland is an important endocrine gland in sheep, and an increasing number of studies have shown that the thyroid gland plays an important role in animal reproduction, but the molecular mechanisms of the thyroid gland in sheep reproduction are poorly understood. In this study, RNA-seq was used to detect transcriptome expression patterns in the thyroid gland between the follicular phase (FP) and luteal phase (LP) in FecB BB (MM) and FecB ++ (ww) small-tail Han (STH) sheep, respectively, and to identify differentially expressed miRNAs (DEMs) associated with reproduction. Bioinformatic analysis of the target genes of these DEMs revealed that they can be enriched in multiple GO terms associated with the reproductive process in animals and in the KEGG signaling pathway. The miRNA-mRNA coexpression network revealed that oar-miR-133 and oar-miR-370-3p may play an important role in sheep reproduction. The results of the dual-luciferase reporter assay suggest a possible targeting relationship between novel-51 and TARBP2. These results provided a novel resource for elucidating regulatory mechanisms underlying STH sheep prolificacy.

Comparative proteomics of ovaries elucidated the potential targets related to ovine prolificacy.

Small Tail Han (STH) sheep, a unique Chinese breed, is recognized for its early maturity, year-round estrus, and prolificacy. However, the molecular mechanism of its high prolificacy has not been fully elucidated. The Proteomics approach is feasible and effective to reveal the proteins involved in the complex physiological processes of any organism. Given this, we performed the protein expression profiling of ovarian tissues during the luteal phase using polytocous STH sheep (litter size ≥2, three consecutive lambings) and monotocous STH sheep (litter size =1, three consecutive lambings) (PL vs. ML), and the follicular phase using polytocous STH sheep (litter size ≥2, three consecutive lambings) and monotocous STH sheep (litter size =1, three consecutive lambings) (PF vs. MF), respectively. Parallel Reaction Monitoring (PRM) was conducted to validate the differentially abundant proteins (DAPs). The tandem mass tag (TMT) quantitative proteomic results showed that a total of 5,237 proteins were identified, of which 49 and 44 showed differential abundance in the PL vs. ML and PF vs. MF groups, respectively. Enrichments analyses indicated that the DAPs including TIA1 cytotoxic granule-associated RNA-binding protein-like 1 (TIAL1), nicotinamide phosphoribosyltransferase (NAMPT), and cellular retinoic acid-binding protein 1 (CRABP1) were enriched at the luteal phase, while TIAL1, inhibin beta-a-subunit (A2ICA4), and W5PG55 were enriched at the follicular phase, potentially mediating reproductive processes in polytocous ewes. Furthermore, six DAPs were verified using PRM, confirming the accuracy of the TMT data acquired in this study. Together, our work expanded the database of indigenous sheep breeds and provided new ovarian candidate molecular targets, which will help in the study of the genetic mechanisms of ovine prolificacy.

Comparative Transcriptomics Identify Key Pituitary Circular RNAs That Participate in Sheep (Ovis aries) Reproduction.

CircRNAs have been found to play key roles in many biological processes and have diverse biological functions. There have been studies on circRNAs in sheep pituitary, and some important circRNAs have been found. But there are still few studies on circRNAs in sheep pituitary with different fecundity. In this study, we obtained the circRNAs expression profiles in the pituitary of FecB ++ genotype Small Tail Han sheep with different fecundity and estrous phases. A total of 34,878 circRNAs were identified in 12 pituitary samples, 300 differentially expressed circRNAs (DE circRNAs) (down: 104; up: 196) were identified in polytocous sheep in the follicular phase (PF) and monotocous sheep in the follicular phase (MF) (PF vs. MF), and 347 DE circRNAs (down: 162; up: 185) were identified in polytocous sheep in the luteal phase (PL) and monotocous sheep in the luteal phase (ML) (PL vs. ML). Cortisol synthesis and secretion pathway (follicular phase) and estrogen signaling pathway (luteal phase) were obtained by functional enrichment analysis of circRNAs source genes. Competing endogenous RNA (ceRNA) network analysis of key DE circRNAs revealed that oar-circ-0022776 (source gene ITPR2, follicular phase) targeted oar-miR-432, oar-circ-0009003 (source gene ITPR1, luteal phase) and oar-circ-0003113 (source gene PLCB1, luteal phase) targeted oar-miR-370-3p. We also explored the coding ability of DE circRNAs. In conclusion, our study shows that changes in the pituitary circRNAs may be related to the response of the pituitary to steroid hormones and regulate the reproductive process of sheep by affecting the pituitary function. Results of this study provide some new information for understanding the functions of circRNAs and the fecundity of FecB ++ genotype sheep.

Effect of melatonin on ATG2B-mediated autophagy regulation in sheep granulosa cells with different FecB genotypes.

Melatonin (MLT) protects cells by reducing reactive oxygen species (ROS) levels, which are key for inducing cellular autophagy. The aim of this study was to investigate the molecular mechanisms underlying MLT regulation of autophagy in granulosa cells (GCs) with BMPR-1B homozygous (FecB BB) and wild type (FecB ++) mutations. GCs collected from small-tailed Han sheep with different FecB genotypes were typed using a TaqMan probe assay, and autophagy levels were found to be significantly higher in GCs with FecB BB than the levels in those with FecB ++. Autophagy-related 2 homolog B (ATG2B) was associated with cell autophagy and was highly expressed in GCs with the FecB BB genotype in small-tailed Han sheep. Overexpression of ATG2B in the GCs of sheep with both FecB genotypes promoted GC autophagy, and the contrary was observed after the inhibition of ATG2B expression. Subsequently, treatment of GCs with different genotypes of FecB and MLT revealed a significant decrease in cellular autophagy and an increase in ATG2B expression. Addition of MLT to GCs with inhibited ATG2B expression revealed that MLT could protect GCs by decreasing ROS levels, especially in GCs with FecB ++ genotype. In conclusion, this study determined that autophagy levels were significantly higher in sheep GCs with FecB BB genotype than the levels in those with FecB ++ genotype, which may have contributed to the difference in lambing numbers between the two FecB genotypes. Autophagy was regulated by ATG2B and was able to protect GCs by reducing the high levels of ROS produced following inhibition of ATG2B through the addition of MLT in vitro.

Progress on the effect of FecB mutation on BMPR1B activity and BMP/SMAD pathway in sheep.

BMPR1B is the first major gene of litter size identified in sheep. However, the molecular mechanism of the FecB mutation that increases the ovulation rate in sheep is still unclear. In recent years, it has been demonstrated that BMPR1B activity is regulated by the small molecule repressor protein FKBP1A, which acts as a key activity switch of the BMPR1B in the BMP/SMAD pathway. The FecB mutation is located close to the binding site of FKBP1A and BMPR1B. In this review, we summarize the structure of BMPR1B and FKBP1A proteins, and clarify the spatial interactive domains of the two proteins with respect to the location of the FecB mutation. Then the relationship between the FecB mutation and the degree of affinity of the two proteins are predicted. Finally, the hypothesis that FecB mutation causes change of activity in BMP/SMAD pathway by affecting the intensity of the interactions between BMPR1B and FKBP1A is proposed. This hypothesis provides a new clue to investigate the molecular mechanism of FecB mutation affecting ovulation rate and litter size in sheep.

The transcription factor RUNX1 affects the maturation of porcine oocytes via the BMP15/TGF-ß signaling pathway.

Bone morphogenetic protein 15 (BMP15) is specifically expressed in oocytes in pigs at all stages from early stages to ovulation and has an important role in oocyte maturation. However, there are few reports on the molecular mechanisms by which BMP15 affects oocyte maturation. In this study, we identified the core promoter region of BMP15 using a dual luciferase activity assay and successfully predicted the DNA binding motif of the transcription factor RUNX1. The effect of BMP15 and RUNX1 on oocyte maturation was examined using the first polar body extrusion rate, a reactive oxygen species (ROS) assay and total glutathione (GSH) content at three time points of 12, 24 and 48 h of in vitro culture of porcine isolated oocytes. Subsequently, the effect of the transcription factor RUNX1 on the TGF-ß signaling pathway (BMPR1B and ALK5) was further verified using RT-qPCR and Western blotting. We found that the overexpression of BMP15 significantly increased the first polar body extrusion rate (P < 0.01) and total glutathione content of oocytes cultured in vitro for 24 h and decreased reactive oxygen levels (P < 0.01), whereas interference with BMP15 decreased the first polar body extrusion rate (P < 0.01), increased reactive oxygen levels in oocytes cultured in vitro for 24 h (P < 0.01), and decreased glutathione content (P < 0.01). The dual luciferase activity assay and online software prediction showed that RUNX1 is a potential transcription factor binding to the core promoter region (-1203/-1423 bp) of BMP15. Overexpression of RUNX1 significantly increased the expression of BMP15 and oocyte maturation rate, while inhibition of RUNX1 decreased the expression of BMP15 and the oocyte maturation rate. Moreover, the expression of BMPR1B and ALK5 in the TGF-ß signaling pathway increased significantly after overexpression of RUNX1, whereas their expression decreased after inhibition of RUNX1. Overall, our results suggest that the transcription factor RUNX1 positively regulates the expression of BMP15 and influences oocyte maturation through the TGF-ß signaling pathway. This study provides a theoretical basis for further complementing the BMP15/TGF-ß signaling pathway to regulate mammalian oocyte maturation.

Epigenetic Regulation of miR-25 and Lnc107153 on Expression of Seasonal Estrus Key Gene CHGA in Sheep.

Pituitary pars tuberalis (PT) plays an important role as the transmission center in the seasonal reproduction of animals. It helps convert external photoperiod signals into intrinsic seasonal reproduction signals. In sheep PT, specific expression patterns of several genes (including short photoperiod-induced gene CHGA and long photoperiod genes EYA3 and TSHß) under different photoperiods are crucial characteristics during this signal transduction. Recent studies have revealed the role of epigenetics in regulating the expression of seasonal reproductive key genes. Therefore, we explored whether microRNAs and LncRNAs regulated the expressions of the above key genes. Firstly, the expression of miR-25 and CHGA showed a significant negative correlation in sheep PT. Results of the dual luciferase reporter assay and miR-25 overexpression indicated that miR-25 could inhibit the expression of CHGA by specifically binding to its 3'UTR region in pituitary cells. Then, expression negative correlation and dual luciferase reporter analyses were used to screen and identify the candidate LncRNA (Lnc107153) targeted by miR-25. Finally, the results of fluorescence in situ hybridization and Lnc107153 overexpression suggested that Lnc107153 and miR-25 were involved in the epigenetic regulation of CHGA expression. However, the expressions of EYA3 and TSHß were not regulated by miRNAs. These results will provide new insights into the epigenetic regulatory network of key genes in sheep seasonal reproduction.

Proteomic Analysis Identifies Distinct Protein Patterns for High Ovulation in FecB Mutant Small Tail Han Sheep Granulosa Cells.

The Booroola fecundity (FecB) mutation in the bone morphogenetic protein receptor type 1B (BMPR1B) gene increases ovulation in sheep. However, its effect on follicular maturation is not fully understood. Therefore, we collected granulosa cells (GCs) at a critical stage of follicle maturation from nine wild-type (WW), nine heterozygous FecB mutant (WB), and nine homozygous FecB mutant (BB) Small Tail Han sheep. The GCs of three ewes were selected at random from each genotype and consolidated into a single group, yielding a total of nine groups (three groups per genotype) for proteomic analysis. The tandem mass tag technique was utilized to ascertain the specific proteins linked to multiple ovulation in the various FecB genotypes. Using a general linear model, we identified 199 proteins significantly affected by the FecB mutation with the LIMMA package (p < 0.05). The differential abundance of proteins was enriched in pathways related to cholesterol metabolism, carbohydrate metabolism, amino acid biosynthesis, and glutathione metabolism. These pathways are involved in important processes for GC-regulated 'conservation' of oocyte maturation. Further, the sparse partial least-squares discriminant analysis and the Fuzzy-C-mean clustering method were combined to estimate weights and cluster differential abundance proteins according to ovulation to screen important ovulation-related proteins. Among them, ZP2 and ZP3 were found to be enriched in the cellular component catalog term "egg coat", as well as some apolipoproteins, such as APOA1, APOA2, and APOA4, enriched in several Gene Ontology terms related to cholesterol metabolism and lipoprotein transport. A higher abundance of these essential proteins for oocyte maturation was observed in BB and WB genotypes compared with WW ewes. These proteins had a high weight in the model for discriminating sheep with different FecB genotypes. These findings provide new insight that the FecB mutant in GCs improves nutrient metabolism, leading to better oocyte maturation by altering the abundance of important proteins (ZP2, ZP3, and APOA1) in favor of increased ovulation or better oocyte quality.

Transcriptome Analysis Reveals Differentially Expressed circRNAs Associated with Fecundity in Small-Tail Han Sheep Thyroid with Different FecB Genotypes.

Litter size is an economically important trait in sheep, and it is a complex trait controlled by multiple genes in multiple organs. Among them, the regulation of lamb number trait by the thyroid gland is a very important part. However, the molecular mechanisms of the thyroid gland in sheep reproduction remain unclear. Here, RNA-seq was used to detect transcriptome expression patterns in the thyroid gland between follicular phase (FP) and luteal phase (LP) in FecB BB (MM) and FecB ++ (ww) STH sheep, respectively, and to identify differentially expressed circRNAs (DECs) associated with reproduction. Bioinformatic analysis of the source genes of these DECs revealed that they can be enriched in multiple signaling pathways involved in the reproductive process of animals. We found that the source genes of these DECs, such as GNAQ, VEGFC, MAPK1, STAT1, and HSD17B7, may play important roles in the reproductive process of animals. To better understand the function of these DECs, we constructed circRNA-miRNA co-expression networks. Dual luciferase reporter assays suggested that a ceRNA regulatory mechanism between circ_0003259-oar-miR-133-TXLNA and circ_0012128-oar-miR-370-3p-FGFR1 may hold. All of these DEC expression profiles in the thyroid gland provide a novel resource for elucidating the regulatory mechanisms underlying STH sheep prolificacy.

Polymorphism and expression of GLUD1 in relation to reproductive performance in Jining Grey goats.

Understanding the molecular mechanism of mammalian reproduction (puberty and prolificacy) will play a part in improving animal reproductive performance. GLUD1 (glutamate dehydrogenase 1) is important for mammalian reproduction, as shown in previous studies; however, its roles in puberty and prolificacy have rarely been reported. In this study, we designed seven pairs of primers (P1 to P7) for cloning and sequencing genomic DNA of Jining Grey goats and Liaoning Cashmere goats. Primer 8 (P8) was designed to detect single nucleotide polymorphism (SNP) of the GLUD1 in both sexually precocious and high-fecundity breeds (Jining Grey, Nanjiang Brown and Matou goats) and sexually late-maturing and low-fecundity breeds (Liaoning Cashmere, Inner Mongolia Cashmere and Taihang goats) by PCR-RFLP (restriction fragment length polymorphism). The real-time quantitative polymerase chain reaction (RT-qPCR) technique was used to detect the expression of GLUD1 in a variety of tissues. The results showed that the A197C mutation was only found in the amplification product of P6. For this SNP locus, only two genotypes (AA and AC) were detected in Nanjiang Brown goats, while three genotypes (AA, AC and CC) were detected in the other five breeds. In Jining Grey goats, the frequency of genotypes AA, AC and CC was 0.69, 0.26 and 0.05, respectively. In Jining Grey goats, AA genotype had 0.54 (P<0.05) and 0.3 (P<0.05) more kids than the CC and AC genotype, respectively, and no significant difference (P>0.05) was found in kidding number between the AC and CC genotype. GLUD1 was expressed in five tissues of different developmental stages. The expression level of GLUD1 in the hypothalamus was higher than that in the other four tissues except during puberty of Liaoning Cashmere goats. In puberty in goats, GLUD1 expression was significantly higher in ovaries than that in the juvenile period (P<0.01). RT-qPCR results showed that the expression of GLUD1 in ovaries may relate to the puberty of goats. The present study preliminarily indicated that there might be an association between the 197 locus of GLUD1 and sexual precocity in goats, and allele A of GLUD1 was a potential DNA marker for improving kidding number in Jining Grey goats.

Characteristics of piRNAs and their comparative profiling in testes of sheep with different fertility.

As a novel class of small RNAs, piRNAs are highly expressed in the animal gonads and their main known role is to inhibit transposon activity for ensuring the correctness and integrity of genome. In order to explore the characteristics of piRNAs in sheep testis and their possible regulatory roles on male reproduction, deep sequencing technology was used to sequence small RNAs and identify piRNAs in testes of sheep. The length of piRNAs in sheep testes showed a unimodal distribution between 26 and 31 nt, with a peak at 29 nt. These piRNAs exhibited obvious ping-pong signature and strand specificity. In the genome, they were mainly aligned to CDS, intron, repetitive sequence regions and unannotated regions. Furthermore, in transposon analysis, piRNAs were aligned predominantly to LINE, SINE, and LTR types of retrotransposon in sheep testes, and the piRNAs derived from each type showed obvious ping-pong signature. The piRNA clusters identified in sheep testes were mainly distributed on chromosomes 3, 7, 15, 17, 18 and 20. The results combining semen determination with pathway enrichment analysis implied that differentially expressed piRNAs between the testes of rams with different fertility might participate in spermatogenesis by regulating multiple pathways closely related to stabilization of blood-testis barrier and renewal and differentiation of spermatogonial stem cell. Taken together, the study provided new insights into the characteristics, origin and expression patterns of piRNAs in sheep testes tissue, which would help us better understand the role of piRNAs in sheep reproduction.

Integrative Proteomics and Transcriptomics Profiles of the Oviduct Reveal the Prolificacy-Related Candidate Biomarkers of Goats (Capra hircus) in Estrous Periods.

The oviduct is a dynamic reproductive organ for mammalian reproduction and is required for gamete storage, maturation, fertilization, and early embryonic development, and it directly affects fecundity. However, the molecular regulation of prolificacy occurring in estrous periods remain poorly understood. This study aims to gain a better understanding of the genes involved in regulating goat fecundity in the proteome and transcriptome levels of the oviducts. Twenty female Yunshang black goats (between 2 and 3 years old, weight 52.22 ± 0.43 kg) were divided into high- and low-fecundity groups in the follicular (FH and FL, five individuals per group) and luteal (LH and LL, five individuals per group) phases, respectively. The DIA-based high-resolution mass spectrometry (MS) method was used to quantify proteins in twenty oviducts. A total of 5409 proteins were quantified, and Weighted gene co-expression network analysis (WGCNA) determined that the tan module was highly associated with the high-fecundity trait in the luteal phase, and identified NUP107, ANXA11, COX2, AKP13, and ITF140 as hub proteins. Subsequently, 98 and 167 differentially abundant proteins (DAPs) were identified in the FH vs. FL and LH vs. LL comparison groups, respectively. Parallel reaction monitoring (PRM) was used to validate the results of the proteomics data, and the hub proteins were analyzed with Western blot (WB). In addition, biological adhesion and transporter activity processes were associated with oviductal function, and several proteins that play roles in oviductal communication with gametes or embryos were identified, including CAMSAP3, ITGAM, SYVN1, EMG1, ND5, RING1, CBS, PES1, ELP3, SEC24C, SPP1, and HSPA8. Correlation analysis of proteomics and transcriptomic revealed that the DAPs and differentially expressed genes (DEGs) are commonly involved in the metabolic processes at the follicular phase; they may prepare the oviductal microenvironment for gamete reception; and the MAP kinase activity, estrogen receptor binding, and angiotensin receptor binding terms were enriched in the luteal phase, which may be actively involved in reproductive processes. By generating the proteome data of the oviduct at two critical phases and integrating transcriptome analysis, we uncovered novel aspects of oviductal gene regulation of fecundity and provided a reference for other mammals.

Characterization of circular RNA profiles of oviduct reveal the potential mechanism in prolificacy trait of goat in the estrus cycle.

The mammalian oviduct is functionally highly diverse during the estrus cycle. It provides a suitable milieu for oocyte maturation, sperm capacitation, fertilization, early embryo development and transportation. While there have been many studies of molecular mechanisms on the kidding number of goats, a systematic analysis by which the underlying circular RNAs (circRNAs) changes in the oviduct related to prolificacy traits is lacking. Herein, we present a comprehensive circRNA atlas of the oviduct among high- and low-fecundity goats in the follicular phase (FH vs. FL), luteal phase (LH vs. LL), and estrus cycle (FH vs. LH; FL vs. LL) to unravel their potential regulatory mechanisms in improving kidding number. We generated RNA sequencing data, and identified 4,078 circRNAs from twenty sampled Yunshang black goats. Many of these circRNAs are exon-derived and differentially expressed between each comparison group. Subsequently, eight differentially expressed (DE) circRNAs were validated by RT‒qPCR, which was consistent with the RNA-seq data. GO and KEGG enrichment analyses suggested that numerous host genes of DE circRNAs were involved in the hormone secretion, gamete production, fertilization, and embryo development processes. The competing endogenous RNA (ceRNA) interaction network analysis revealed that 2,673 circRNA-miRNA-mRNA axes (including 15 DE circRNAs, 14 miRNAs, and 1,699 mRNAs) were formed, and several target genes derived from the ceRNA network were associated with oviduct functions and reproduction, including SMAD1, BMPR1B, IGF1, REV1, and BMP2K. Furthermore, miR-15a-5p, miR-181b-5p, miR-23b-5p, miR-204-3p, and miR-145-5p might play important roles in reproduction. Finally, a novel circRNA, circIQCG, was identified as potentially involved in embryo development. Overall, our study provides a resource of circRNAs to understand the oviductal function and its connection to prolificacy trait of goats in the differentiation estrus cycle.

Expression characteristics of piRNAs in ovine luteal phase and follicular phase ovaries.

PIWI-interacting RNAs (piRNAs), as a novel class of small non-coding RNAs that have been shown to be indispensable in germline integrity and stem cell development. However, the expressed characteristics and regulatory roles of piRNAs during different reproductive phases of animals remain unknown. In this study, we investigated the piRNAs expression profiles in ovaries of sheep during the luteal phase (LP) and follicular phase (FP) using the Solexa sequencing technique. A total of 85,219 and 1,27,156 piRNAs tags were identified in ovine ovaries across the two phases. Most expressed piRNAs start with uracil. piRNAs with a length of 24 nt or 27-29 nts accounted for the largest proportion. The obvious ping-pong signature appeared in the FP ovary. The piRNA clusters in the sheep ovary were unevenly distributed on the chromosomes, with high density on Chr 3 and 1. For genome distribution, piRNAs in sheep ovary were mainly derived from intron, CDS, and repeat sequence regions. Compared to the LP ovary, a greater number of expressed piRNA clusters were detected in the FP ovary. Simultaneously, we identified 271 differentially expressed (DE) piRNAs between LP and FP ovaries, with 96 piRNAs upregulated and 175 piRNAs downregulated, respectively. Functional enrichment analysis (GO and KEGG) indicated that their target genes were enriched in reproduction-related pathways including oocyte meiosis, PI3K-Akt, Wnt, and TGF-ß signaling pathways. Together, our results highlighted the sequence and expression characteristics of the piRNAs in the sheep ovary, which will help us understand the roles of piRNAs in the ovine estrus cycle.